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^ ■ Abstract 

a^ 

l— ■ A new approach to quantum gravity is presented based on a nonlinear quantization scheme for 

\^ . canonical field theories with an implicitly defined Hamiltonian. The constant mean curvature foliation 

^D I is employed to eliminate the momentum constraints in canonical general relativity. It is, however, 

^I ' argued that the Hamiltonian constraint may be advantageously retained in the reduced classical sys- 

■^ ■ tem to be quantized. This permits the Hamiltonian constraint equation to be consistently turned 

qh, into an expectation value equation on quantization that describes the scale factor on each spatial 

5" I hypersurface characterized by a constant mean exterior curvature. This expectation value equation 

{^' augments the dynamical quantum evolution of the unconstrained conformal three-geometry with a 

transverse traceless momentum tensor density. The resulting quantum theory is inherently nonlinear. 

Nonetheless, it is unitary and free from a nonlocal and implicit description of the Hamiltonian oper- 



\^ • ator. Finally, by imposing additional homogeneity symmetries, a broad class of Bianchi cosmological 



models are analyzed as nonlinear quantum minisuperspaces in the context of the proposed theory. 
PACS numbers: 04.60.Ds, 04.20.Fy ll.lO.Lm, 04.60.Kz 

1 Introduction 

The term 'quantization' is ill defined. Ideally, one starts with a quantum theory and may wish 
to classicize it. However, until the ultimate quantum theory of everything is found, one continues 
to 'prepare' a classical system for quantization under a consistent scheme. Further justifications 
of the quantized system may then be sought. In the canonical approach, this requires a choice of 
phase variables. Gauge or other auxiliary degrees of freedom, if any, must be factored out from the 
redundant phase space so that only the physically significant degrees of freedom are quantized. The 
corresponding type of variables are often referred to as dynamical variables whose identification from 
a given theory is not always straightforward and can be ambiguous [I|- Thus, a central question to 
ask in quantizing a classical system is: what is the right point of departure for quantization? Once 
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this is decided, the actual quantization may then be addressed in terms of its physical consequence 
and mathematical expediency. The problem of quantization is particularly perplexing for gravity. 
Starting with the Dirac- Arnowitt-Deser-Misner (Dirac-ADM) formalism of general relativity |2j , the 
implementation of canonical quantization has been impeded by the presence of the Hamiltonian and 
momentum constraints that signal the redundancy of the conventional geometrodynamical variables. 
While the momentum constraints are clearly related to the invariance of the formalism under spatial 
diffeomorphisms, a gauge interpretation of the Hamiltonian constraint is more evasive. In particular, 
the Dirac algebra satisfied by these gravitational constraints fails to be a Lie algebra of any gauge 
group. 

Two divergent strategies have been adopted in dealing with the gravitational constraints. The 
first strategy aims to exploit them as first class constraints to generate the quantum dynamics of the 
spatial geometry through the Dirac constraint quantization 12 lEl El EI • In essence, this scheme turns 
both Hamiltonian and momentum constraints into quantum annihilation equations for the physical 
state functional. The second strategy is based on the reduced phase space method, whose basic 
theory and history may be found in |S1 [HI El ^D 13 El ^2 El • When applied to canonical general 
relativity this strategy advocates the elimination of all gravitational constraints in order to obtain a 
nonvanishing effective Hamiltonian that generates the evolution of certain unconstrained geometric 
variables '16[ 1171 1181 I19j . The reduced system would be amenable to the conventional Schrodinger 
quantization in principle. 

The Dirac constraint quantization strategy gains support from the confirmation that the Dirac 
algebra guarantees the hypersurface-independence of the dynamic evolution and hence the general 
covariance of any parametrized field theories |2Ul I21j . Together with the interconnection between 
the Hamiltonian and momentum constraints [22 this property suggests that, under Dirac constraint 
quantization, a hypersurface-independent quantum evolution of gravity could be constructed without 
the need for gauge fixing. As the the Hamiltonian is quadratic in the canonical momenta of the spatial 
metric, this strategy results in the Wheeler-DeWitt equation j231 I24j . which is a functional Klein- 
Gordon equation for three-geometries that has no obvious probabilistic interpretation. While this 
non-unitary quantum evolution may lead to interesting predictions such as the tunnelling phenomena 
in certain models of the early universe |23], considerable efforts, notably by Kuchaf [^HlilZlilHl) have 
been made to find alternative set of gravitational variables in which Dirac quantization yields the 
Tomonaga-Schwinger type equations |291 I3UJ that generate unitary quantum evolution of geometry. 
For a broad class of midisuperspace models with two commuting spatial isometries, the symmetry- 
reduced general relativity can be shown to be equivalent to a parametrized scalar doublet field theory 
in 2-dimensional Minkowski space jHJ. However, it is not clear whether Kuchaf 's variables can be 
found in full general relativity |S21 ■ Further discussion on Kuchaf 's method may be found in ^53^ and 
references therein. 

The strategy of applying the Schrodinger quantization to a reduced gravitational phase space by 
means of gauge fixing was initially explored by Dirac JHIJ and later by ADM P] . The ADM procedure 
amounts to choosing four out of the six spatial metric components as special spacetime coordinates. 
The gravitational constraints are solved (and are thereby eliminated) for the canonical momenta of 
these coordinates, leading to a nonvanishing 'true' Hamiltonian. The remaining two spatial metric 
components with the corresponding momenta are then regarded as dynamical variables. Note that 
this procedure does not by itself single out any preferred dynamical variables. The analogy that 
ADM drew between the gravitational and electromagnetic fields led them to introduce a background 
Minkowski metric, with respect to which the transverse traceless part of the general spatial metric 
components and their momenta are treated as dynamical variables to be quantized. However this 
scheme relies on an ad hoc reference metric at the expense of loosing covariance. 



In this regard, a much more satisfactory gauge fixing scheme that singles out a 'spin-2' part of 
gravitational field was developed by York J35ll36l07j . Built upon previously works by Lichnerowicz, 
Brill and Choquet-Bruhat 38 , 39 , 4QJ . York defined a decomposition of general second rank symmetric 
tensors of any weight into transverse traceless, longitudinal and trace parts, in a conformably covariant 
fashion. The application of this powerful scheme to canonical general relativity is most revealing if 
the time-slicing endows each spatial hypersurface with a constant mean exterior curvature. Under 
this constant mean curvature (CMC) foliation, the spatial metric is written in terms of the scale 
factor (of the volume element) and a conformal metric. The latter is used to decompose the metric 
momentum tensor. The resulting longitudinal part of the metric momentum tensor must vanish to 
satisfy the decoupled momentum constraints. The trace part of this tensor is proportional to the 
mean exterior curvature which labels the spatial hypersurfaces like time. This 'exterior time' may 
also be considered as a manifestation of the geometric carrier of information about time originally 
suggested by Baierlein, Sharp and Wheeler [IJ. Therefore, the transverse traceless part of the metric 
momentum tensor, along with the conformal metric, naturally carries dynamical degrees of freedom 
and can be interpreted as gravitational waves jSZl- The scale factor is to be determined by solving 
the reduced Hamiltonian constraint as an elliptic equation, i.e. the Lichnerowicz equation |42l I43j . 
Furthermore, it has been shown by Isenberg, O'Murchadha, York and Nester that York's conformal 
treatment is capable of accommodating scalar, electromagnetic, spinor fields [44_ and indeed 'virtually 
all' known physically relevant forms of matter field j45j . 

The clear geometrical meaning and physical interpretation of York's identification of the dynamical 
variables of general relativity raised the hope for quantum gravity by quantizing the conformal three- 
geometry with transverse traceless momentum. The diffeomorphisms of the three-geometry and the 
longitudinal part of the metric momentum (nullified by solving the momentum constraints using the 
CMC condition) may be regarded as gauge variables. In contrast, no simplicity seems to exist in 
dealing with the remaining Hamiltonian constraint. If this constraint is also solved for the scale 
factor before quantization as per the ADM procedure, then a nonvanishing Hamiltonian for canonical 
general relativity may be constructed, which turns out to be the volume integral of the universe. As a 
nontrivial elliptic equation must be solved, this Hamiltonian is a nonlocal functional of the dynamical 
variables and extrinsic time and is only defined implicitly. This somewhat disturbing feature has 
been noted by Choquet-Bruhat, York, Fischer and Moncrief j43 ( I46| IT7| I4SJ . It presents at least a 
technical obstacle to quantization as issues such as factor ordering ambiguities are hard to tackle with 
an implicit and nonlocal Hamiltonian |49j . 

Under this state of the affair, alternative quantization schemes should not be ruled out where 
York's identification of the dynamic part of gravity may be utilized without having to solve the 
Hamiltonian constraint prior to quantization. The purpose of this paper is to explore one such 
possibility. Specifically it is proposed that the scale factor may be treated on an equal footing with 
the lapse function as a Lagrangian multiplier as befits their roles in formulating the kinematics of 
the spatial hypersurfaces. These non-dynamical variables are conceptually akin to Dirac's surface 
variables in his analysis of the generic surface kinematics of the parametrized classical field theory 

The rest of this paper is devoted to developing the above new strategy for canonical quantum 
gravity, with emphasis on its physical interpretation and geometrical basis. In section |21 the previously 
proposed nonlinear generalization of a finite-dimensional physical system based on a quantum action 
principle j52l I53j is extended to a field-theoretical description. This generalization is motivated by 
the 'implicit nature' of the gravitational Hamiltonian on the CMC foliation. Before proceeding, it is 
necessary to lay out the essential prerequisites for classical canonical gravity. This is given in the first 
half of section |31 with a focus on the elimination of the momentum constraints leading to a reduced 



system of the pattern discussed in section |2j The second half of section El addresses the quantization of 
the 'CMC-reduced' canonical general relativity in sectional using the nonlinear quantization method 
developed in section |2l The corresponding quantum action generates the quantum evolution of the 
conformal three-geometry augmented with two elliptic equations that determine the scale factor and 
lapse function. In order to establish the essential framework of the proposed theory with a view 
to more complete treatments, some of the discussion given in this paper remains formal. As in the 
development of other quantum gravity ideas, a 'toy model' in an even simpler mathematical setting 
will help to grasp the salient features of the new approach. Of course, one should always bear in mind 
that the usefulness of this way of extracting guidance is limited [511 ESI- Despite this caveat, it is 
demonstrated in section 0] that the Class A Bianchi cosmologies and the Kantowski-Sachs universe 
can be treated as nonlinear quantum minisuperspaces of the proposed theory. The conclusion with 
discussion on future work is given in sectional Below units oi c = h= IGvrG = 1 are adopted. 

2 Nonlinear quantization of a field theory with an 
'implicit Hamiltonian' 

Let M be a three-dimensional space with coordinates x* {i = 1, 2, 3). With respect to a preferred time 
r consider the canonical evolution of a field (j) and its conjugate momentum vr, each with n components, 
namely (j)"" = 4>"'{x,t) and tTq = 7ra(a;,r) (a = 1,2,- --n). The time-dependent Hamiltonian density 
h = /i((/), (/)^j, • • • ,7r,7r^j,--- ,r) is assumed to be explicitly and locally defined. Here and henceforth, 
the subscript comma and over dot denote differentiation with respect to a spatial coordinate and time 
r respectively. The canonical field equations are generated by varying the action integral 

S[(P,tt]= ff iTT-<j)-h\d^xdT (1) 

where ir ■ (p := vTa 0", with respect to vr and (j) to be 

in terms of the Hamiltonian 

H= / hd^x. (3) 

Jm 

The treatment of a canonical field theory using an implicit Hamiltonian can be developed as 

follows. In terms of an arbitrary time t an action equivalent to (^ is 

S[(p, vr, a, N] = [ [ {tt ■ dfCp - a dtT - NH] d^x dt (4) 

J Jm 

where 

n = -a + h (5) 

and T = r(t) together with two Lagrangian multipliers a = cr{x, t) and A^ = N{x, t). Note that 7i has 
no explicit dependence on the parameter time t. The integral 

S= /" ad^x. (6) 
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is introduced for later use. 



Now put t = T but retain r and A'^. The action Q) becomes 

S[(l), TT,a,N]= 11 i'K-^-a- Nn\ d^x dr. (7) 

This form suggests a generahzation from the restricted form of W in @ a generic 7i = 
TC{(j),(p^i,- ■ ■ ,Tr,TT^i,--- ,a,a^i,--- ,t). For the purpose of the subsequent discussion, it suffices to 
assume that this expression is also exphcitly and locally defined. With this generic TC, the variations 
of the action ((Tj) in vr, (j), N, a generate, respectively, the equations 

■ 6H 6H 

W = (9) 

l.f.O (10) 

where 

H = f NTid^x. (11) 

Jm 

Equations Q govern the evolution of (j) and vr. Through H these equations are coupled with © and 
(|1U|) that can be regarded as algebraic or elliptic equations determining a and N respectively. For a 
to be interrelated as the 'energy density' of the system, it should be nonnegative. Clearly, the above 
equations reduce to to ^ for the 'simple case' where © holds, with © and (|Tn|) yielding a = h 
and A^ = 1. Given a generic 7i, although in principle @ may be solved for the 'true' nonvanishing 
Hamiltonian density a, this will bring about a high degree of nonlocality into the reduced canonical 
theory. Such an implicit and nonlocal Hamiltonian may in turn cause severe conceptual and technical 
problems for quantization. 

The standard canonical quantization of the system with the action (^ yields the Schrodinger 
equation 

i^ = H * (12) 

for the state functionals ^ = ^[i;^] in terms of the Hamiltonian operator 



H= /" hd^x. (13) 

Jm 



where h is obtained by substituting vr — > vr = — i^ into h with a suitable choice of factor ordering. 
The Schrodinger equation (|12j) can be shown to follow from the quantum action integral |531 I56j 

5q[^,**] = /"{Re(*,i^) - (*,H*)} dr (14) 

using the inner product of any state functionals ^i and ^2 in terms of the functional integration 

(^i,*2>:= I ■^\-^2'D<P (15) 



in an appropriate sense. The operator h and other physical observables are assumed to be Hermitian 
with respect to this inner product. The expectation value of any Hermitian operator O with respect 
to ^ is defined, in a standard manner, as: 

(O) := ^ff^. (16) 

By the same token as the steps leading from (^ to (0]) , an action equivalent to ()14() can be written, 
using the Lagrangian multipliers a{x,t) and N{x,t), as 

5Q[*,^*,cj,iV]= /"{Re(^,i^)-S(*,*)-(*,i^^)}dr (17) 



where 



with 



H= / NTid^x (18) 




n = -a + h. (19) 

Here Ti. is the operator of TC given in ©. For a generic TC = TC{(p, 4>^i-,- ■ ■ , vr, vr^j, • • • , vr, cr, ct j, • • • ) 
as discussed above that admits a corresponding Hermitian operator ?Y, the use of this operator in the 
quantum action H17() will, under variations with respect to \1', its conjugate, A^ and o", generate the 
following generalized form of quantum evolution equations 

i^ = H^ (20) 

{ft) = (21) 

= (22) 

respectively, up to a (physically irrelevant) overall time dependent phase of ^ . 

The Schrodinger type equation H2U|) thus describes the unitary quantum evolution of ^[(j)]- The 
'Hamiltonian operator' H depends on a and N ^ for which H21|) and ()22() . must be solve. As these two 
equations are nonlinear in the state functional, the overall quantum evolution is nonlinear. Although 
the introduction of nonlinearity into quantum gravity is relatively new, the nonlinear modification 
of quantum mechanics has long been receiving considerable attention in addressing the measurement 
problem |571 1581 IHU] . Interestingly, the energy incurred by the nonlinear term in certain 'dynamical 
reduction models' can be shown to be comparable to the gravitational self-energy of the measurement 
apparatus [^Zl- This leads one to suspect that gravity might be related to quantum evolution and 
measurement in a nonlinear manner. For further discussion on nonlinear quantum theories, see e.g. 
|6Ul 1611 I62j . In view of this, the nonlinear quantum formalism presented in this section has been 
developed with a view to quantizing general relativity. It will be demonstrated in the next section 
that the action integral for classical canonical gravity does indeed reduce to the form of (O under the 
CMC foliation. The nonlinear quantization scheme may therefore be employed for quantum gravity 
as discussed below. 



3 Canonical general relativity on the CMC foliation and 
its nonlinear quantization 

In this section the essentials of the Dirac-ADM formulation of canonical general relativity are first 
recapitulated. The theory is then gauge-fixed by choosing the spatial hypersurfaces to possess a 
constant mean exterior curvature. The resulting 'CMC-reduced' formulation is readily organized 
into the form of the generalized canonical field theory discussed in the forgoing section that may be 
quantized in a nonlinear fashion. 

Under the ADM 3+1 split of spacetime, each spatial hypersurface arbitrarily labelled by time t has 
a spatial metric g with components gij in coordinates x* {i = 1, 2, 3). The inverse (i.e. contravariant) 
metric g~^ has components g^^ . In terms of the lapse function N and shift vector X with components 
X^, the spacetime line element takes the form 

ds^ = -N^ dt^ + gij (dx* + X' dt){dx^ + X^ dt). (23) 

Here indices are lowered or raised using gij and g"^^ respectively. 
The extrinsic curvature tensor K has components given by 

^'' = ^ ^-9iJ + ^-^J + ^r,i) (24) 

where the semicolon denotes a covariant derivative using the Levi-Civita connection of g. This connec- 
tion gives rise to the intrinsic scalar curvature TZ. Of particular importance for the present discussion 
is the mean exterior curvature defined as the trace of K by 

/C := g'^K.j (25) 

In terms of K and /C the 'metric momentum' tensor density p is given by 

p'^ = -fi {K'^ - g'^IC) (26) 

where fi = (det (7)2 is the scale factor for the spatial volume element. 

The well-known Dirac-ADM action integral for canonical gravity takes the form 

S[g, p,N,X]= ft {p-g-NH-X -J} d^x dt (27) 

J Jm 

where p ■ g = p^-'gij, X • J = X^Ji in terms of the Hamiltonian and momentum constraints 7i and Ji 
given by 

W = - fi'^{gik9ji + 9u9jk - 9ij9ki) p'^P^^ - ^/U (28) 

Ji = -2gij p^'^-k = -2 gij p''',k - 2 gij,k P^'" + 9jk,i p'^ (29) 

respectively. 

The roles of N and AT* as Lagrangian multipliers result in the Hamiltonian constraint equation 

H = (30) 

and the momentum constraint equations 

Ji = 0. (31) 



York's treatment of canonical gravity begins by expressing the metric g as 

9 = ^S (32) 

in terms of a confornial factor 93 and a conformal metric 7 satisfying a 'normalization' condition. This 
may be done by fixing its volume scale factor ^^ = (det 7) 2 [121150] or its scalar curvature IZ.^ whereby 
regarding 7 as a Yamabe metric [031 EH HE] ■ E-g- 



^^ = 1 (33) 

or 

7^^ = ±l,0 (34) 

depending on the topology of the spatial hypersurface M. The present discussion does not depend 
explicitly on the choice of such normalization conditions. However, boundary terms in various inte- 
grations over M are dropped at will for the sake of mathematical convenience. That is to say, either 
a compact M or an adequate fall-off of the relevant fields over M at infinity is assumed. 

An interesting result of York's is that the momentum constraints decouples from the Hamiltonian 
constraint under the CMC foliation. Furthermore, the momentum constraint equations (|31|) are 
satisfied if and only if the longitudinal part of the metric momentum p vanishes, leading to the 

expression 

2 

p = ip^'^w + -ip^]C^^^^.^f (35) 

o 

where ro is a transverse traceless tensor density with respect to 7 (HHl EH • 

The Hamiltonian constraint equation gives rise to the Lichnerowicz equation |38[ 1361 1461 I5UJ 

3 11 

A^(/? - ^ rV^ - g n^V + g liklji ^'^^^^ /^7 ^ ^'^ = (36) 

regarded as the elliptic equation for the conformal factor ip if H30|) is also to be solved. The present 
analysis will however defer this process by retaining 7i in the action. Since 

M = V^^/^7 (37) 

in the subsequent discussion the conformal factor p will be eliminated using 

_i 1 

(P = H-y 6 /U 6 . (38) 

Following |46l I47j one writes 



p-g = ^-j + -ICfi (39) 

where zu-'j = w^^ 'jij. This, together with the fact that (|3H) has been solved by (|35)) and an integration 
by parts, reduces the Dirac-ADM action ()27() to the form 

S[-f,vj,^i,IC,N]= ff lvj-^--fifC-Nn\d^xdt. (40) 



Here the CMC-reduced Hamiltonian constraint Tl = Ti{'y,zu, fi,T) is obtained by substituting (|32j) 
and H35() into ()28() and has the form 

n = -l nT^ + iM-^^ikiji w'^w''^ - nil (41) 

o 

Q 

= -g/iT"2 + /^"Sifc7j7ro*-'cCT''''-^7^73At3 +8//^fe^6 A^(/i^"6/i6) (42) 



where A^ is the Laplace-Beltrami operator on scalar functions with respect to 7 (following the sign 
convention of [131 EOl as opposed to that of [HZl |^)- Note that the action (|in)) has a structure 
analogous to that of (HJ). To bring (|4n|) to the form of 0, where a, 4> and n may be identified as fj,, 
7 and w respectively, the time parametrization t = t := o/C is adopted, yielding the action 

S[-f,w,n,N]= 11 {w^-iJi-Nn]A-^xdT (43) 



I M 

where 



and 



w^d''x-Vol-H\dT (44) 

J 

H= f NHd^x (45) 

Jm 



Vol= / /ud'^x (46) 

JM 



representing the 'volume of the universe'. It follows that 



5H _4. ^ _ _r3 
5^ 3 

_ 4 

~ 3 



AiV - iv| ^ t2 + ^-27,,,7^, w'^w^' + i 7^| (47) 

5 _ 5 _i In 

A^l g T^ + f^~'^liklji ro^^w^' + - 7^^/i^i/i-i - - fi^6^-l A^(^f i/x^) i (48) 



where A is the Laplace-Beltrami operator on scalar functions with respect to g. 

In complete analogy with Q, © and (fTU]) . the variations of the the CMC-reduced action H43|) 
with respect to w, 7, A^ and /i generate the canonical type field equations 

. SH 5H 

for the conformal metric 7 and its transverse traceless momentum density zu, supplemented with 

n = (50) 

and 

1 + ^ = (51) 

as the elliptic equations for the scale factor /j, and lapse function N, respectively. 

Proceeding with the nonlinear quantization scheme formulated in section |21 one now considers 
a state functional ^ = ^t[i] depending on the conformal metric 7 modulo diffeomorphisms |49j . 
The inner product of any two such state functionals ^1 and ^2 are given in terms of an appropriate 
functional integration of the form 

(*i,*2):= /^*i^2^7 (52) 




over the equivalent classes of the conformal metric 7 related also by difFeomorphisms. Corresponding 
to the classical action ()44() . one further writes down the quantum action for the above ^t[7]) f^ix,T) 
and N{x, r) as 

Sq[^,^*,H,N]= fiRe{^,i4')-Vol{^,^)-{^,H^)\dT. (53) 

This action in turn generates the following Schrodinger type equation 

i^ = H^ (54) 

for ^ (up to an overall time dependent phase), coupled with the elliptic type equations 

{n) = (55) 



= (56) 

which serve to determine n and N, respectively. In the above equations, the operator 7Y is obtained 
by substituting w ^ w := —i-^ into ()42() where the term quadratic in vu naturally becomes the 
Laplace-Beltrami operator in the space of conformal three-geometries. It is worth noting that the 
above procedure still carries through if the conformal factor f is chosen in place of //, with H56|) 
slightly modified as ^ becomes an operator in terms of 7. 

The Schrodinger type equation H54|l provides a unitary quantum description of the conformal 
three-geometry. Through the transverse and traceless nature of the corresponding momenta, it is 
gravitational waves that is quantized. The quantum evolution depends on /i and N which are deter- 
mined by the elliptic type equations H55|) and (|56() . These two equations make the overall quantum 
evolution nonlinear and nonlocal in the the state functional. Nevertheless, the Hamiltonian constraint 
operator is locally defined. Its kinetic part is quadratic in the momentum operator and should be 
represented by the Laplacian in the space of conformal three-geometries. 

4 Nonlinear quantum minisuperspaces 

The nonlinear framework for quantum gravity presented in the previous section is complicated by 
the presence of the Laplacian of // and N as well as the field-theoretical nature of the problem. 
This complexity can be avoided in the cosmological approach with finite degrees of freedom. At this 
oversimplified but concrete level, one expects to gain intuition about the full theory for being able to 
carry out explicit analysis. 

For this purpose, the Class A Bianchi type models and Kantowski-Sachs universe |651 I66j will be 
treated as nonlinear quantum minisuperspaces of the full theory under discussion. These cosmological 
models are described by the lapse function N(t), length scale factor R{t) as well as two additional 
functions /9+(t) and f3-{t) to allow for the dynamical anisotropy of the spacial hypersurface at any 
parameter time t. The spacetime line element takes the following common form 

ds^ = -N^dt^ + R^{e^% e'e^ (57) 

for i,j = 1,2,3, where e* denote some basis 1-forms of the spatial hypersurface and /? is a traceless 
matrix with elements given by 

f3ij = diag(/3+ + V3/?_, /3+ - ^/3/3_, -2/3+). (58) 
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The basis l-forms e* satisfy the structure equation 



de* = ^ C)k e^ A e*^ (59) 



with certain model-specific structure constants C^jk |65j . Evidently, the spatial metric g may be 
written as 

9 = R^l (60) 

where 

1^3 = {e^% (61) 

satisfying the normalization condition (jSHJ- Hence fi = B? . 
It follows from (gl, (HSJ, (EHl) and dnSJ that 

with the over dot indicating differentiation with respect to t, and 
?3 



w^^ = — diag 



Therefore 



(/3+ + ^/3/3_)e-2^+-2^'^-, (/?+- V3/3_)e-2^++2V3^-, -2/3+ e^^+ . (63) 



tu • 7 = ^ (/?+ + /?-) = ^+/9+ + ^-/3- (64) 



where 

^± = ^/3±- (65) 

With the choice of time t = t and the substitution (as simple volume normalization) J^j-d^x — > 1 
the action (jlH)) reduces to 

5[/3±,ro±,/x,Af] = /"|ti7+/3+ + ti7_/3_-/x-iVH|dr (66) 

where 



82"^ 



2 



1 



^ = "8/^-^ + ^-^^^. (67) 

with ro^ := vd\ + tuL 

As no factor ordering ambiguity arises, the cosmological model may now be quantized using the 
standard substitutions: w± — > ■dj± = — i9/3_,_ and ro^ — > w'^ = —{dp + d'i_). These operator are to 
act on a wavefunction ^ = ^(/3±, r). The inner product of two such wavefunctions ^1, ^2 is simply 

00 /•oo 



(^i,*2>:= / / ^^^2d/3+d/3_. (68) 

As per preceding discussions, the quantum action of (|H^ can be written as 

Sq[^,^*,H,N]= /|Re(^,i^)-^(*,^)-7V(^,?i:^)}dr. (69) 



This generates the Schrodinger type equation 



1 



i^ = ^^|^-^^7^^!^M/ (70) 
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for ^(/9±,r) (up to an overall time dependent phase), and following two equations 

-^.r^ + |^-.M^.)^0 (71) 



l-A^^^^ + ^ + ^/^-^W^=0 (72) 



(ct^) 1 2 

2V ^3 

algebraic in /^(t) and N(t). Equation ()7U() describes the unitary quantum evolution of the anisotropy 
(i.e. 'truncated gravitational waves') of the cosniological model in a mechanical fashion. In particular, 
it involves a 'Hamiltonian operator' consisting of the sum of a kinetic term and a potential term. 
However, the presence of the supplementary conditions (|71|) and (|72|) introduces nonlinearity as well 
as nonlocality into the description. In this manner, equations H7U|) - (|72() constitute a consistent system 
of nonlinear integro-partial differential equations. Such equations can be solved at least numerically 
in a fashion similar to [SI] where a computational approach to the nonlinearly quantized Friedmann 
universe is detailed. In the case of the Bianchi I model, the properties TZ.y = and (tz7^) being a 
constant enable one to find exact solutions to H7U|) - H72() using the analytical method developed in |53j . 

5 Concluding remarks 

A discussion has been given of a nonlinear quantization scheme for canonical field theories with an 
implicitly defined Hamiltonian. Two Lagrangian multipliers are involved in the formulation: one 
representing the value of the energy density and the other effecting this value as an implicit function 
of the canonical variables allowing for explicit time-dependence. This structure has been shown to 
naturally arise from canonical general relativity on the constant mean curvature foliation, providing 
a basis for a nonlinear quantum theory of gravity. As the scale factor (or conformal factor) and 
lapse function remain unquantized, they act as Lagrangian multipliers in a quantum action principle. 
In particular, the scale factor represents the effective positive energy density. The true dynamical 
degrees of freedom for general relativity, as identified by York, are carried by the conformal three- 
geometries with transverse traceless momenta and are quantized nonlinearly in the proposed theory. 
York's exterior time, namely the constant mean exterior curvature, or physically the (local) rate of 
expansion of the universe, is also adopted here as the preferred time. This work has set the stage 
for more explicit and detailed investigations. Indeed, work is underway to transcribe the present 
geometrodynamical description to a 'connectodynamical' description so as to assimilate powerful 
functional techniques offered by the loop quantum gravity approach. As the presented framework 
features inherent nonlinearity as well as nonlocality in the state functional, it is of interest to explore 
the resulting physical consequences. This will form a subject for future research, together with other 
issues such as the inclusion of matter interaction. On this last point, it is reasonable to anticipate 
that much of the essential methodology of the present work for pure gravity to be carried over, owing 
to the versatility of the constant mean curvature analysis in accommodating matter fields. 
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